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SUMMARY AND INTRODUCTION 


This report presents the results of experimentation with a single 
jet-augmented flapped airfoil for generating gusts of high reduced 
frequency for the me a 2 e = 

‘OD 
in a wind tunnel. The experimental data from a single jet flap in 
conjunction with the theoretical assumption of superposition, are 
used to postulate combinations of trailing vortices from an assumed 
array of muitiple jet flaps for the production of an oscillating 
cross-sectionally-uniform gust. The objective of this report is to 
study and understand the effects of different parameters such as flap 
deflection angle, flap frequency, dynamic pressure, and the position 
at which the induced velocity is measured, on the uniformity of the 
generated turbulence across the wind tunnel test section and finally 


to combine these system variables to produce sufficiently uniform gusts. 
Three areas of work were the main points of concentration: 


1.) Experimental testing in both the static and dynamic cases 
with a single jet flap or vortex generator in the wind 
re Cigi@lsale 

2.) Superimposing the results of the experimental data into an 
array of four generators in an attempt to produce a uniform 
cross-sectional gust. 

3.) Comparison of the gust model designed using the experimental 
data and superposition with the model using a theoretical 


analysis of a single trailing vortex and superposition. 


ial al 


The testing was accomplished in the Princeton University 4' X 5' 
wind tunnel with jet augmentation of a simple plain flapped airfoil. 
The jet augmentation increased greatly the magnitude of induced velocity 
generated, and the use of a flapped airfoii decreased the size and inertia 


of the moving components. 


The experimental results consist of measured values of induced 
velocity due to a single half span vortex generator. Both static and 
dynamic data were taken. In the static case, the magnitude of induced 
velocity was dependent on certain variables such as flap deflection angle 
sensor vane position in the test section and free stream velocity. 
Statically, many combinations of parameters were observed and recorded. 
Thus, the partial derivative of induced velocity with respect to each 
parameter was separately observed. For the dynamic data, induced 
velocity was dependent on-Strouhal Number. Bi-directional as well as uni- 
directional gusts were generated and were compared. An analysis of the 
effects of high frequency cancellation and vortex filament movement 
illustrated a similarity in peak to peak magnitude in both types of 


generated turbulence. 


The second area of study was the use of the experimental data in the 
hypothesis of an array of generators that would produce the desired uniform 
gust model. Superposition was used to pattern the symmetrical array of 


trailing vortex filaments into this desired model. 


Finally, a similar array was developed from a trailing vortex 
filament in an inviscid flow field. This theoretical model was compared 
to the model using the experimental data and superposition previously 


designed. They were in close agreement. 


sir 


LIST OF SYMBOLS 


As | Slot Area; As=on 5 (Chee) 
AR Aspect Ratio of jet augmented airfoil ; 
b | Span of the vortex generator (ft) 
G@ Average chord of gust generator (ft) 
oA Drag component of the sensor vanes 
Cd. same as oG 
O% 
Con Minimum drag coefficient 
Cs Chord length of the flap (ft) 
a Total lift coefficient of the vortex generator 
Cn Blowing coefficient of the vortex generator; 


Lo ay, 
7 “9.9 / 


Dr De Drag on the fore and aft airfoils of the sensors, 
respectively (1b) 


DYNAMIC Time variation with flap deflection. 


a Perpendicular distance of a point in space to a vortex 


‘core which is inducing a flow at that point (ft) 


|, - Width of the slot (achk2s) 
As Streamline deflection distance in 2-D smoke tunnel 
testing - (inches) 
ii Moment of inertia or the vane potentiometer system (slug-ft ) 
Lk’ Reduced frequency - nondimensional frequency of flap 
oscillation 


. Kal 
K rod 
Eafe 


pe LAL 


De Lh 


Full scale reduced frequency 
Model reduced frequency 


SistinicaninGiscance Of TUllsalreraLt, used in computing 


reduced frequency (ft) 
Induced velocity perpendicular to the generator chord ( direction) 


Moment arms of the fore and aft airfoils of the sensors, 


respectively (1b) 


Significant distance of model used in computing reduced 


frequency (ft) 

Mass flow rate of jet air through the blowing slot; 
wh; = e AY; 

Point in test section to be observed 

Pressure at upper boundary of the jet sheet (lb/ft”) 

Pressure at lower boundary of the jet sneet Gb/tta) 

Dynamic pressure of the test section (1b/ft”) 

Radius of curvative of the na sheet (ft) 


Reynold’s Number; 7 Vok 
Jie 


Area of the vortex generator airfoil Gan 


Constant maximum flap position. 


Area of the airfoils on the sensing vanes (ft) 


Jet ejection velocity; always considered sonic in the 


present study (ft/sec) 

Induced velocity parallel to the generator span (~ direction) 
Free stream wind tunnel velocity (Fé/s Lc) 

Induced velocity at point P (ft/sec) 


Magnitude of (ft/sec) 


Streamwise distance from the leading edge of the test section (ft) 


=< 


Distance from the left wall of the test section looking 


downstream (feet) = - 


Nj 


Vertical distance from the base of the test section looking 


downstream (feet) 


Keo Effective angle of attack (deg) - angle measured from the 


ZeLOuisibmian Male 
| Induced angle of attack at the vortex generator (degrees) 
ena pe Angle at which the shed vortex core will trail off behind the 
flap (Degrees 
Oo Induced angle in the wake at infinity (degrees) 
oCre Stall angle of attack (degrees) 
OG 2-D circulation at upper boundary of jet sheet (ft /sec) 
ae 2-D circulation at lower boundary of jet sheet (ft/sec) 
Lf 2-D circulation caused by emitted jet sheet (tt /sec) 
Pr 


Circulation about the shed vortex core (ft/sec) 


r. Total circulation causing pressure lift (ft/sec) 

ie Total circulation causing momentum lift (ft? /sec) 

eS Derivative holding all other variables constant 

i. Magnitude of flap deflection either positively or negatively 


(degrees ) 


‘Se Sensor position in streamwise (X) direction as well as «+ 


in the cross-sectional plane (Y-Z) 
E oo Same as (so 


5 Non-dimensional damping coefficient of sensor vanes 


av anlat 


a) Angle of incidence of the relative wind to the airfoil 


surfaces of the vane sensors (radians) 


ORF Angle involved in the determination of the influence of a 


semi-infinite vortex of a point in space (degrees) 
YA Wave length of the full scale gust to be simulated (ft) 
A mop Wave length of the model gust (ft) 
e Density of air; constant (1b/ft ~) 


om Torque produced by the aerodynamic loading on the airfoils 


of the sensors (ft-1b) 
GJA/c Full scale gust frequency (hz) 
C.) 5 Flap frequency of oscillation (hz) 
Cipro Model gust frequency (hz) 


Con Natural frequency of vane sensors (rad/sec) 


gs bales t 


DESCRIPTION 


Gust simulation in the 4 X 5 foot Princeton University wind tunnel 
was conducted using a blown flap gust generator. No attempt was made to 
improve the design of the gust generator based on experimental observations. 
However, changes were made in the sensors during the trial runs to improve 
the data output and to increase experimental accuracy. The individual parts 
will be discussed separately including the changes in the sensor GOnstrucL ion 


fons this trial run period. 
Wind Tunnel 


The tunnel test section (Figure 1), nominally four feet by five feet, 
is forty-two inches high and fifty-nine inches wide. This Princeton 
facility has the least tunnel turbulence of all the wind tunnels at the 
Forrestal Campus. Its dynamic pressure is variable over a range of tunnel 
velocities up to approximately 100 mph. The maximum tunnel free stream 
velocity used in the experiments was 100 ft/sec which is below the upper 
limit where tunnel turbulence tends to be amplified. The test section 
has a streamwise length of 80 inches, thus allowing surveys se three positions, 


two, six and ten chord lengths downstream from the jet flap model. 


Inherent in any wind tunnel is its turbulence caused by wall surface 
irregularities and friction (shear) stresses. This turbulence has its owm 
characteristics of frequency and amplitude which must be filtered out of 
the final data in an attempt at some degree of accuracy. The frequency 
of the turbulence was recorded in the initial trial runs as 0-50 lz 
which was well above the testing range of the experiment and was easily 
distinguishable and filtered. The magnitude of the turbulence varied 
between one and two degrees of free stream deflection in a steady state 


condition, and was filtered out of the data presented. 


a 


JET FLAP VORTEX GENERATION 


The jet flap was designed specifically to increase the lift 
Sapacrlaty of a particular airfoil section. Since the (; is increased, 
the circulation also increases and hence the strength of the trailing 


vortex far exceeds that of a simply flapped airfoil. 


To create a Sinusoidal gust the airfoil section was chosen as a 
NACA 0012 symmetric airfoil. The flap, (Figure 2) was taken to be a 
5c flap for reasons discussed in APPENDIX EZ. Blowing over both sides 


of the flap was necessary to insure a symmetrical gust profile. 


This trailing vortex generator was a hal? span model mounted to 
the tunnel floor plate to facilitate oscillating the flap. Mechanically, 
the flap is oscillated by means of a connecting rod attached to a driven 
fly wheel as in Figure 3. The wheel has different radial attachments 
for the connecting rod. By increasing the radial position the actual 


flap deflection can be increased both positively and negatively. 


The drive wheel is powered by an electric motor through a 10 to l 
ratio reduction gear box with a right angle power transmission. The 
motor was rated for maximum power output for 27 volts with 69 amperes of 
current. Throughout the entire range of tests the voltage never increased 
over 27 volts with 10 amperes; therefore, the maximum power output was 
never approached. All powered components are controlled from a central 


console as shown in Figure }. 


Blowing over tne flap wes-accomplished by the introduction of va 
015 in,widtn slot over the full twelve inch semi-span cf the airfoil. A 
large plenum chamber was designed into the flap construction to feed air 
to this slot. The Schramm Compressor was able to produce a plenum chamber 
pressure of approximately 30 psi, enough to choke the slot and give spanwise 
uniformity to the jet sheet. Line loss was almost 50 psi from the cutput 


Or the compressor to the point of induction into the plenum chamber. 
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FIGURE |. 


DOWNSTREAM VIEW OF TEST SECTION. 
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VORTEX GENERATOR. 


JET AUGMENTED FLAP, 


FIGURE 2. 
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ELECTRIC’ D.C. MOTOR, GEAR BOX, 


DRIVING WHEEL, AND CONNECTING ARM. 


DRIVE MECHANISM FOR JET FLAP; 


FIGURE 3. 
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FIGURE 4. 


CONTROL CONSOLE AND WIND TUNNEL. 
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FIGURE 5. 


INS TRUMEN TA TION 


Air stream following vanes mounted on low inertia potenticmeters 
in Figure 5 were used tO measure the steady state and time varying 
response of the free stream air to the sinusoidal oscillation of the jet 
flap airfoil upstream. The dynamics of the vanes and their design 
characteristics are discussed in Appendix IV. Adequate speed of response 
and no overshoot in time response are two desirable characteristics of the 
vanes. After many different types of vanes were designed and tested, the 
vanes shown in Figure 5 were used. They have Low aspect ratio conical 
bodies with laterally-protruding wings at leading and trailing edges. 
By inducing an impulse perpendicular to the free stream in a steady 
state condition it was observed that with an aft moment arm of 3 inches 
and a forward arm of 1.5 inches, the vanes are critically damped at all 
tunnel velocities. This is advantageous for accuracy in data recording. 
The natural frequency was in the vicinity of 30-35 Hz which is weil beyond 


the range of experimentally-induced frequencies. 


The probe shows the configuration of the potentiometers 
and vanes in the free stream. One potentiometer is mounted laterally 
placing its vane in position to record any vertical component of 
deflection in the stream. The vertically mounted pot and vane will record 
the horizontal conponent of disturbance. By vectorially adding these 
components we arrive at the velocity vector induced at that point. If we 
assume that this velocity vector is perpendicular to a radius to the 
vortex core then this represents a good method for experimentally locating 


the vortex core position laterally in the tunnel. 


Electrically excited by a 15 volts source the potentiometers were 
attached to a Sanborn heated stylus recorder and a Brush pen recorder. 
With maximum frequency capability of 60 eps and 100 cps, the Sanbern and 
Brush recorders, respectively, could easily accommodate even the high 
tunnel turbulence frequencies. Calibration was accomplished by mechanically 


varying the vanes through a known angular distance. Brush and Sanborn 





recorders had 1 mm displacement on the recorders equal to 1.5 degrees 


of vane deflection for the former and 1.3 degrees for the latter. 


EXPERIMENTAL PROCEDURES 
The experimental testing was accomplished over a period of six weeks. 
A number of assumptions had to be made to simplify reduction of the data. 
Incompressible flow and constant temperature were assumed. With a 
number of variable parameters, the most efficient means of surveying 


cross sectional areas of the tunnel was devised. 


The variable parameters consist of the following: 


as 


{ 


Tunnel dynamic pressure which alone varied the yoru. 


Cie ie lea pede mileceina om, 


1.) 

es) op 

3.) Frequency of flap deflection, i. 

4.) Position at which the induced velocity is measured in 
test section, Jy Eve , 


YY 


The technique consisted of the following steps: 


Pipsc ts. Cae telloyvaric: 
\ 
The specific *;- value. 
o 
A specific vane position in the tunnel, sypyé. 
The initial tunnel dynamic pressure, q. 


Sea re eee e vane 
The initial flap frequency, -<. 


The range of these parameters is very important in gust simulation. 
The various combinations can change the magnitude, direction, and span- 


Wise distortion of particular gusts. 


16: 


For the present discussion the following magnitudes of the 


parameters were employed. 


The initial condition consisted of 


7 ¢ 


First, a static case was run for each new «? and vf. The blown flap 

was manually cranked, slowly, through a single cycle. This created 

a good static trace. Then the flap was oscillated through a range of 

fae, > nz, hk hz, 8 hz, and 16 hz, with the run time of 10 to 15 

seconds. Next the Vowas lowered to 75 ft/sec and the static and dynamic 

cases were run again. Finally, the ‘> was fixed at 50 ft/sec and the process 
was repeated. The ve and t¢ changes were made for all Seen a specific cross 
meecion Of tunnel at x = 2.0c downstream from ne model. After this cross 
section was surveyed, the vanes were moved to the x = 6.0c downstream 

Section and finally, the x = 10.0c downstream section. The YE magnitude 

was then changed to accommodate a variation of + 45° . After the entire 
tunnel was surveyed again at Ne = + 5° : the Sr-was changed to + 60 
for the third and final run through the tunnel. During the Bea = 30 
Bievey, approximately eleven to thirteen evenly spaced points were surveyed 
miseach cross section. For the Se =+h5° and + 60° , only four of the 
centrally located points in each section were observed. This reduction in 
survey points does not decrease the accuracy of the experiments. Induced 
velocity magnitude variations are still discernible and comparisons with 


different parameter regimes can still be made. 


A few runs were made without blowing augmentation for comparison with 
the augmented generator. And finally, data were taken with the connecting 
- 
fimeshortened, to oscillate the flap unidirectionally from Lh a Oe to 


See 60° for dynamic observations-discussed later. 


DISCUSSION OF EXPERIMENTAL RESULTS 
STATIC: 


In the static case the experimental trends were correlated with the 
theory. In Figure 8 a typical central position of Dig was chosen and 
observed. Diffusion of the vortices downstream is negligible from a 
theoretical as well as experimental point of view. Westmore (Reference 29) 


relates a variation of (Swith downstream distance as follows: 


By substitution it is observed that, with the tunnel "q's’ in the present 
discussion and the short length of streamwise distance, diffusion down- 
stream is negligible. Even with the additional turbulence introduced by 
the tunnel walls, no dissipation can be observed over this distance. [In 
faeu, 2S is observable in Figure 7, there is a strengthening in las the 
probe is translated downstream. This is most likely due to the decrease 
in 4 and increase in cos ©). for positions farther downstream when Seeore 


(Figure 19). 


Having assumed that all of the vorticity is wrapped up in a single 


trailing vortex, the profile of downwash (Reference 18) is as follows: 


a 


Although this graph is induced angle of attack with distance X, 1t gives 

a good indication cfr the actual variation with distance of the vorticity 
as a combination of circulation and jet sheet components to the total 
memo.) Th1S, therefore, 1S a good represencation of the strength of the 
trailing vortices. Here then we can see that as X increases it is obvious 
that at some central position [ov p{ will be greater than at an even farther 
ayaié downstream. Observing the actual tendencies in Figure 7 we see 
that this hypothesis holds true in all cases. In the x = 2.0c position, 
the large value of and small cos Gp magnitude, due to core position, 


) 
at x = 6.0c downstream. This discussion represents the static case with 


will overcome this initial |’ at x = e+*hence, the value of iJp will peak 


the flap deflected positively, (toward the centerline of the tunnel). It 
-also introduces the interference tendency far downstream due to wall 


effects that are not noticeable near the flap. 


When the flap is directed in a minus direction, (toward the closer 
tunnel wall), h will increase with x as will cos £). 5 dais atone 
effects tend to cancel in the generation of lps Again wild ie) as eae 
Figure 6. Theoretically, the general trend for [2¥.| should be decreasing 


in value with downstream distance. This is normally true (Figure 19). 


In the set of Figures 8 to 16, the entire cross-section is viewed. 
Both magnitude and direction of the ‘J: vectors make the position of the 
core in the static condition easily located. The solid vectors are for 
as dp and the hatched vectors for Sp Pee eon Guidemanicn directa On 
differences between the hatched and solid vectors illustrate the movement 


of the equivalent vortex core position across the tunnel. 


a 
In any cross-section, Dyno movement obviously decreases SR 


with increasing \ proving experimentally that 


OS p T= ~ 


13 


Velocity ana ':( changes alone can be observed by comparing Figures 8, 
9,and 10. These are typical arrays of cross-sectional induced velocity 


vectors with varying free stream qelocye (te. Thus, as is’ evident from 


J We = 
3 Ver 


« a ® e \ ° 
The general trend is toward an increased iS for an increase in rer 


however, the partial derivative value is not strong. Hence, this variation 


these figures, 


is somewhat erratic. Experimentally, (Figures 17 and 18), this solution 


agrees with the theoretical evaluation in Appendix I. 


a 
Flap deflection changes, OJ, have the most obvious affects on 


Of all the parameters. As can be seen in Appendix f[, 


a) 1 oF 


All three add a positive increment to the partial derivative such that 
9 Ws 
ere 

SS 


A comparison, in progression, of Figures 8, 11 and 14 will give a typical 
= 


— 
—_ 
a 


view of the change in magnitude of the induced velocity as a function of of 
Figure 20 gives a typical graphical picture of this progression. It can 
be Observed that the 1133] (Se) seems to be of the order of an increasing 


mower function. It appears then that in fact the partials 


are es 
O 2 fe } a é have 
‘a 


greater influence at higher values of Of 
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FIGURE 7. COMPARISON OF THE MAGNITUDE OF INDUCED 
VELOCITY AT A CENTRAL POSITION AT VARYING 
DOWNSTREAM POSITIONS WITH THE THEORETICAL 
INDUCED VELOCITY AT X= IN THE WAKE OF 
THE JET FLAP. 
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FIGURE 8. COMPARISON OF CROSS-SECTION SURVEYS OF INDUCED 
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DISCUSSION OF EXPHRIMENTAL RESULTS 
DYNAMIC: 


The dynamic experimental data for the majority of cases agree with the 
theoretical trends discussed in Appendix I. In all regimes, the decrease in 


induced velocity amplitude with high frequency is evident. 


Induced angle of attack holding all the variables but W-, constant was 
recorded as a function of reduced frequency (Strouhal Number), K. As can 
be seen from Figures 22 and 23, magnitude peaks in Wp arise at intermediate 
values of K for those sectional positions investigated. As K is increased 


these data tend to show the same decreasing trend as Figure el. 


The rapid decay of magnitude of induced velocity as a function of K in 
Figures 21 - 23 at high frequency shows strikingly the dynamic interference 
affect of high frequency oscillation of the airstream. Some of the high 
frequency data, however, does show a reduction trend about a biased value 


of (Sp; this is due to the vortex filament movement downstream. 


In the augmented case, as is stated in the theory of Appendix I, at 
high K, the magnitude of iJ, , for the sinusoidal symmetrically-oscillating 


flap does level off to a constant value in the range SS Zz up 
He ile AX &F 
i 


foerigures 21 = 23, 2 comparison was made of the induced flow 
generated without blowing augmentation. It is obvious that the trends Oi 
peaking with intermediate values of K and decay with high values of K 
Seer here, too. But as before, the use of blowing creates a great increase 


in magnitude even’‘at high K values. 
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Experimental results (Figures 2 : 27) indicate unidirectional flap 
oscillation decreases the cancellation of the harmonic components which 
appear in the symmetrical case. However, the theory proposes that this 
high frequency cancellation should be exactly the same”as that for the 
symmetrical deflection. The theory and experiment disagree at this point 


Since the aS] does not diminish symmetrically around the mean value. 


Vortex filament motion as described in Appendix I will only change 
the value of the bias such that at (S-= OO the | Sarhen the fa lamers 
assumed motionless,is less than [u5>] when the filament is in motion. 
This fact is still of less importance than the previous disagreement 


between the calculated and the experimental cancellation. 


Although further experimentation at higher K is necessary to verify 
the experimental decay of (Oo with K, it would seem impossible to create ~ 
gust by means of an oscillating jet flap and trailing vortex system for 


otvrouhal Numbers much in excess of 2.5. 
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NUMERICAL EXAMPLE 


To compare the proposed theory of Appendix I with the actual ex- 
perimental data, a typical set of parameters is chosen to be analyzed 
in the static regime. | 


The following are assumed: 


a 
(1) dyes » Such that 


X = 6.0c downstream 
Y = 30" 
Z = 30" 
oe O | 

(2) dQ¢ = 45° , static in the positive direction 

(3) Wee = 100 ft/sec with a choked 
a) 

ng can S . 8 
From Table I: 
ie 
(1) vc, = 3.49 
(Ce eae 


(3) Filament angle (90° 3) Swaner? 
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From the Biot-Savart Law 
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The value of h when calculated from Figures A and B yields 
h = 1.88 ft 


Also assuming this value of h, then 


5 BONS i 
SIA) SON a (2) 
Thus, the cos Ip = .923 

Substituting these values into (1) and assuming (J behind the flap 


along the jet sheet is Fe aR NOC oe induced velocity at point P is 


as follows: 


rh 


f : }+ (234 4). 





— () 


Be) i > WEE C Hees 8 


h (3) 


and finally 


This value represents a single vortex in free air. Now the experimental 
wind tunnel results can be seen in Figure il, where the magnitude of 


induced velocity is 


Wp a MOno Ft/seEc 


Soe 
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ins COrrelat.on Holds true wi Chri phe range experimental error, for 
the great majority of surveyed points. This would lead one to believe that 
some parameters in the wind tunnel are possibly changing the wall effects 
and cancelling them. In Heyson's report (Reference 12) the affects of both 
Spen and closed wind tunnel tests sections are analyzed. The overall effects 
of these two are opposite in sign. It also makes reference to the possibility 
Bemearcelling any wall erfects by introducing a slotted or porous tunnel 
Section. Something in the tunnel, which is seemingly obscure at this time, 
is simulating this porous tunnel effect. Further study is necessary to 


discover the exact problems. 


ho 


PROPOSED MODET, 


The ideal gust disturbance would be a two-dimensional gust. It would 
seem reasonable that such a disturbance pond be achieved by the use of 
superposition of vortices. This report proposes just such a system. Although 
previous research has generated different methods such as the Texas A&M 
pulsing jet for longitudinal gusts (Reference 27) and Giannini Controls 
Corporation oscillating wall jet gust simulator (Reference 6), the 
oscillating airfoil method used at NASA Langley (Reference 7) for generating 
@ sinusoidally oscillating airstream seems best for a two-dimensional gust 


generator. 


As Reference / indicates, many problems arose in generating high 
frequency gusts in an attempt to simulate high speed flight conditions. 
Appendix I demonstrates the inertial load problems that would indicate 
selection of a simple flap mechanism with blowing augmentation to produce 


Similar, if not greater, gust magnitude capavilities. 


Purely vertical or horizontal gusts can thus be generated by a 
superposition model. By oscillating the jet augmented-flaps either in 
phase or out of phase by neo » these gusts can be produced. The vortex 
array should be patterned to generate as near a constant cross-sectional 
Piiobeas possible. Thus, the gust seen by the model being tested will be 
nearly constant in a spanwise direction. As will be seen later, with the 
flaps oscillating out of phase by 180° , the most uniform gust will remain 
on the centerline while varying in magnitude and direction. However, when 
oscillated in phase, the filament movement will cause an unwanted siift in 


this most uniform gust position. The complex velocity potential is, 


ejp--t Lb ea (2) 


men 


a1) 


mor & Single vortex, where 


Se ee 
Z+Xty ee 


a = location of the vortex in the plane. 


eance ie (2) = - ( Lind at a=o0O 
A 


by substitution of (2) into (3) we arrive at 


or 


Now 
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By comparing (5) witn (6) we see that 
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(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 
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Substituting we arrive at the radial velocity 


Vy- = © (11) 
and the tangential velocity 
a 
ie ae (12) 
Oo” amr 


This equation illustrates the inverse proportionality between the radial 


distance and magnitude of the vortex-induced velocity. 


Since the proposed system has an array of four vortices its velocity 


potential is assumed to be the following, 


‘ i ; , cnet | : { : x . \/ 
fs (2)- ~eE j,[ 2 -(y, thy) ss ide LJ? 68 “a ) | 
| a > s ac ae eee oy \ 
tee by [z r (x mn 4) ite le AE 7 (x- oD) . (13) 


The array would appear as in Figure 28. 
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puch an array would produce a purely vertical oscillating gust. 

In producing a gust perpendicular to the flap chord, the dynamic 

movement of the vortex filaments at positions farther downstream 

will all be out of phase. However, even in a high frequency mode 

the lateral velocity components will all cancel in the central area 

of the tunnel. The position of the most uniform spanwise gust velocity 
will not move as a function flap deflection. This presents the best use of 
Secillating jet flap position, frequency and phase of oscillation, to date. 
If the jet flaps were mounted on the side walls of the tunnel a pure 


horizontal gust could be generated. 


From the complex potential previously described can be derived the 
Bonjugate velocities in the x and y directions. 
? an © Omi 
OF: 


Where u = velocity in x direction 


v = velocity in y direction 


By Se Sorc uso liets (eke Sas exists 
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Using these equations, the theoretical magnitude and direction of the 
induced velocity can be shown. The movement of the vortex filaments is 
involved in this analysis. The equivalent vortex positions at x = 6.0c 
were assumed. This provides a means of comparison for this case with 


the actual experimental value. 


The dimensions of Figure 28 will be assumed since the experimental 
data will be of that configuration. With the gust generation devices 
mounted to the floor and roof of Due sono Unters Laps ware Symchroni2 ed 
180° out of phase with the bordering flaps as shown in Figure 28. This 
will provide as even a vertical gust as possible. Figure 29 illustrates 
theoretically the magnitude and distortion of :Swhich must be tolerated 
when testing under such a configuration. This distortion, for a 2-foot 


span model, is approximately 40% from mid-span to tip. 


Figure 30, is the comparable experimental result due to super- 
position. Here again, spanwise distortion, in vertical w is important. 


It is approximately 40% over the same 2-foot span model. 


dy 
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Creating lateral gusts would seem at first, easily accomplished 
by changing the phase of the flap oscillations. Figure 32, is an 
illustration of this configuration. The generators in the array are 
all of the same phase. However, the filaments all moving in phase 


back and forth are ignored. 


In Figure 32, experimental superposition shows that the center of 
the array of filaments as appears in Figure eal will move a certain 
distance to the left as the filament particles move downstream. Hence, 
at maximum static deflection, at x = 6.0c downstream, a very erratic 
pattern appears in the centerline position. This is intolerable when a 
relatively constant cross-sectional gust is required. In the production 
of these gusts, both lateral and vertical, combinations of parameters 
could possibly be found and programmed to induce a uniform gust flow 
Svyer the entire span of the model. In addition, a half span model could 
be mounted to the tunnel floor with the previous vertical turbulence, 


Peerepy Simulating a lateral gust. This,however, is left for future work. 
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CONCLUSIONS 


There are some basic conclusions that can be drawn from the 


discussion and analysis of the experimental data and theoretical 


calculationsas presented previously. is 


1) 


2) 


ai) 


Since only the flap and not the entire airfoil is being 
oscillated, the inertial loads at high frequencies can 


be decreased considerably. 


By properly programming the vortex generators, the array 
of airfoils can produce a uniform cross-sectional gust 


‘Woe une Vertical or lateral directions. 


In the dynamic case the flap was oscillated both bi- 
directionally and uni-directionally. Theoretically, the 

former should produce an oscillation in induced velocity 
symmetrically centered about zero. At very high frequencies 

the dynamic components tend to cancel themselves. Experimentally, 
This was Substantiated. Theoretically, the latter should be 


centered about a finite value such that when (Jo = 
si 


Ce 
(Sy > 4 ve 


vey ly jie ae 


However Che experimental results dn this Series of tests did nev 
fully agree. The experimental dynamic peak to peak value did not 
decrease symmetrically about the constant value. However, 
Peete lesmcOnmmcs pond LOouougouual Numbers dm-excess of 1.5 
must be observed and more data taken to substantiate this 


deviation from the theory. 


MY) 
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pi 


Vortex filament movement is a factor in determining the value 
of the constant value of Wp when wr) ¢ a ct 

With the present apparatus, the steady state value for the bi- 
directional oscillation will be displaced from the zero 
position. This, however, does not influence the fact that the 
time varying oscillation is cancelled as if the core position 
did not vary. Core movement in the uni-directional oscillation 
again simply biases the steady state value, but again, had no 
effect on the dynamic cancellation. The experimental data 
confirms these theoretical observations. Hence, in the 
production of a very high frequency, sinusoidal, time varying 
gusts, the vortex generator is of little use. However, the 
apparatus can be operable for a low to medium frequency band- 
width while taking into account the rapid magnitude decay at 


the upper band frequencies. 


Attempts to account for the presence of the Panne Dereael by 

image vortex systems indicate that little or no wall effects are 
necessarily accounted for in the calculation of ob, In the 
numerical example this fact was shown for an average test point 
where it can be seen that the test section actually simulates 
Cpe Lousesceulen, and, Loeoretically, needs Jittie or mo tunnel 
wall correction. Further investigation is needed to discover the 


reasons for this trend. 
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APPENDIX I 
JET FLAP ANALYSIS 


Important aerodynamic parameters of the jet flap were fixed due to 
the problem itself. First, an airfoil of a symmetrical cross section 
was desirable to accommodate the changes in circulation throughout a =< 
tO tog range which would in turn produce symmetrical gusts. The NACA 0012 
section was chosen for conveniencé because templates with a nine inch 


chord were already in existence, thus simplifying construction. 


Once the chord was fixed, the flap size had to be determined. 
Preliminary two dimensional work,described in Appendix II, was done in 
the Princeton 2" X 36" smoke tunnel to arrive at a possible ©, for the 
3-D model. Knowlton, in Reference 14, discussed methods of determining 
m6. for a 2-D airfoil section oy measuring the streamline deflection 


of an upstream position. 


Various tunnel dynamic pressures, flap deflections, and Cu. Sere 
observed with varying flap chord ratio. The tests were conducted at 
30-50 ft/sec. Thus, the data were for the range of Reynold’s Numbers 
ee < Me Ze and, hence, this regime corresponds very well 


with the 3-D model. 
VORTEX FILAMENT 
From 3-D elementary airfoil theory without blowing argumentation, 
LL oa ce q S ai ) 


Seesubstituting for a constant chord model, 


o { 
aes | L \/ a ise ce 
Le a ee { J f 2) fy C (T-ia) 


From the simple circulation theory of lift and the Kutta Joukowski 


theorem, i | 7 
ee (1-2) 


such that lift per unit span is equal to the preduct cf alr density, 


free stream velocity, and circulation. 


By equating (la) and (2) we arrive at 


a ee wiki 
C, (4 io") be = bp Vee | (1-3) 


Eliminating terms by dividing through we arrive at a value for 


Cimeulabion as a: function of ©, Veo andc. 


) 
Hence, from (3) 


j' ~ Cee 


=) 
a 


(r-4) 
Since 


Fo a. fy 
C2 (x, Sr] ak a ~ O (1-5) 


shal that G is a function of angle of attack and flap deflection angle 
With no blowing augmentation, we realize that when the angle of attack 


is zero then lift coefficient depends only on flap deflection or 
a CS a ae 
© = 9 ere Qk KO (1-5a) 


Her the model used in these experiments, the flapchord is nearly 
half the total chord. Let us assume that the change in angle of zero lift 


, : C ji 
is approximated by half the Ss or the change in angle of zero lift is, 


he: 2 7 
oe = WS one (1-6) 
pn? 
A 
Since, for any non-augmented airfoil,', is a constant, 


ne 


then, im A 
[' i “ts 9s Ve ; (I-7) 
pei im 


Thus, for an unaugmented, plain-flapped airfoil at c<. = 0 
(a 


until the flap deflection reaches a value of twice the stalling angie 


Of attack, or when 


> 45 : | 7 
Sf = ee ee (I-8a) 
The average induced angle of attack, with a fully attached flow 
7 
is dependent only on Ops IRlahbiss 
8 if a 
isthe coi <A O os (I-9) 
Vee 
In the case of the blowm flap augmentation (Reference 19), the lift 
a 


©eefficient is a function of both Jef and oe where 


@ . 4 (5; ) ~ (I-10) 


due to the introduction of Ut into the analytic determination of on 


in Reference 20. 


Hence, circulation for the blown flap is proportional to the product 

a fr 
of ihe and ©, , with the latter a function of both OF ays = LOUS 7 Her 
this case 
| 


es | 
/ — : { Sy, Cu, Ve (I-11) 
Now by definition . 


C. =e Oo i | 


(I-12) 


ai 
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in that the product of the mass flow rate of jet air and its Local slot 
velocity is nondimensionalized by the mass momentum over the wing. Thus, 
Cu mein fact, a lift coefficient in addition to the on or circulation 
lift coefficient. \/7 is constant throughout the present experiments; 
Bence, CC 


= 


promit ReTerencele>,, che inerement in LiftGeas by reaction to the 


is dependent only on ve 


vertical component of the jet momentum and the increases in wing 
Sereulation lift from the introduction of asymmetry to the airfoil by 
the jet and the energy injected into the boundary layer to retard 
separation. Both increases in lift will be assumed to add a component 


of to the trailing vortex shed from the tip of the 3-D wing. 
JET SHEET 


From thin airfoil theory, the jet must first be considered in 


inviscid flow. Figure 33 gives an illustration of the problem. 





The high velocity air jet is issued from the. leading edge of the flap, 
leaving the flap surface at the ‘¢ angle. Both the flow inside the jet 
and surrounding the jet are irrctational. Total head pressure is greater 
in the jet than in the free stream. Figure 34 illustrates a segment of 


the jet downstream. 


De 





7 
os _ \ 
Joa / ‘ a 


u( re 
- Rr, hj 
gece “44 
Since the flow has a discontinuous velocity gradient, we can consider 
the jet surrounded by two opposite signed vortices each proportional 


to the strength of the local flow external to the jet. 


Again from Reference 25, the pressure at the boundaries is 
continuous, i.e. total head is constant across the boundaries. From 


Bernoulli's equations, 
ua 


Ne mer ) 
Pr + y @ UL | a iF One. US Geen) 


yy ee 
ie + a OS = a L® v (Tai 


Since irrotational flow is concerned 


} a ie / 
Ue ie i. y 4) re i (IR : a (I-15) 
| ~ ~ 3 


‘ 


By involving mean velocities in the discussion and rearranging algebraically 


we can arrive at the conclusion that 


ae 
aie —— Ty) @ Vv 
P~pe 7 R : 
a — (I-16) 
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where 


ment . 
V4 p b (uv, +0;) (lai7 


Again this assumes that the jet extends to infinity and thrust is 


totally recovered. y 


The pressure differential described infT-16)is due to the difference 
in the elemental strength of the upper and lower boundary vortices. 


By integrating the two into one vortex by assuming 


: ° 
| = O 
— —. 
a 


un Gene) 


we will arrive at a vortex sheet with a strength equal to the 
difference in strength of the previous two, located at the centerline 


cece jet. Thus, | 
Oj ds re (we ds, = LS ds, (I-19) 


Following Spence's progression in theory, the strength of the vortex 


jet sheet per unit span is 


y (I-20) 


But since L>c from(I-18}and on, must remain constant then \/- must 


J 
approach ©° as 


VY. & ( Ye. (1-21) 
. — oF loo C | 
(Ez 
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and the pressure differential per unit span can be assumed from (I-16) 


to be y 
wma of e 
fe ca 


eed 


R a (1-22) 


Now the total lift is the sum of circulation lift and momentim 
i Laka 


JJ (I-23) 


4 
Where le is the corrected circulation lift for the blowing augmented flap, 


and 
CO 
= 
ie ram a) ¥- as (T-2h) 
J ; 
Since 
: tho 
Y= 
Olas eo Cu (I-25) 
J 2R 
g os 
Then the following thin airfoil approximations are made for small J = 
and «. 


1) the vortex distribution representing the jet is located from 


I< X¥<© along the real axis; 


2) yy (x — Le Nee Cu 4 “x) 


where y (x) is deflection of the jet below the horizontal line; 


(y) ee 762 for Ce | 
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Oe Fig uke os 


vy) = Vu ¥(y) -y V40 Cu yl) 


By following through with integration of(I-25a),taking the function to 


where 


(I-25a) 


teoelimits, and solving for CL 8 total of the force components will 


give values of . 
2.54 Cu © $326 Un 4,196 

S. ae ies ye Va | ieee 
and 


, OP Fie, 406 
a = Vaya eye Curt , Aah r oe GEE 


With further refinements Spence's equations are displayed in McCormick 


(Reference 20) as follows: 


a | eee oa ( |” 
dle [a (113) ye all as (i=28) 


3 &. , 
and | = : | / . 
oo 7 or ¢ yas ly (1-29) 


C/- 
fee 


Reference 20 was used for actual interpolation between the above two 


These equations hold true for wo NOwevchemcie eraphaacal plou am 


equations. 
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For total () the following is true: 


ee: ($\ de sie oh ~ ; (I-30) 


since the present flap chord is in the region of 
O<°SK <1.0 
an interpolation was necessary, and a numerical /, was obtained. 
Table I lists the numerical values for all experimental regimes calculated 
using McCormick's interpolation. Also ene actual value of *%<¢ 
is listed for varying values of fy, and) ee Numerically, here it can 
be seen that (,, (or V2) has little effect on the average downwash angle 


\ 


at infinity; but that it is nearly linearly dependent on ee 


a | 
Thus CO og ash 0 (Tea 
oO VO 


oO 
Ack ea 
ea eee 


with a fairly constant positive value. This in turn shows that the 





(1-32) 


while 


position of the vortex filament dowmstream, as shown in Figure 36, is 
not a function of °..? but only of Sr. Also this leads to the conclusion 


that »~. and cos re Gem nGenencene OL hei 


For analytic purposes assume the value of Ca given in (29) to 





be correct. “A 
Thus , C “ (J GC ‘ 
1K C ———— 
¢ 5) (1-33) 
es 7 c a 
where es = is mid eras ( 6} 








TABLE I 
Vee (Ft [ec >= 50 75 
Cu — 3,33 Lae 
oF 
30° 3.4 2.62 
nS 5.1 3.8 
60° 6.8 oe On 
BO 25.7 26.2 
BS 38.6 37.9 
60° 51.5 50.6 
Bo” Ds) Deon 
45° 7.9 aa 
60° 10.6 10.4 
ASSUMPTIONS: 


aR val 


TER pC 


100 


Bh 
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VORTEX FILAMENT POSITION FOR SINGLE VORTEX 
l0.0c 


5=20° 


5,=0° 
6.0c7 


DOWNSTREAM 


20c-4 


A i nee eel AEE eee, lS eee Gee eee ee aS ee eee 


FLAP TRAILING EDGE 


FIGURE oC. 


Since from (4) 


a ease) 


ae ‘ V. 
p= §;[¢ he tsa | et (1-35) 


where 
Ly: (I-35a) 


From the law of Biot and Savart, the induced velocity anywhere in the 


system due to the circulation in a single trailing vortex is, 


IW, = Iisa dl 


7p 7 (I-36) 


From Dwinnell (Reference 4), in the steady flap condition, the total 


induced flow at a point "P" due to the entire trailing vortex is 


LI > . 6 di Wp (I-37) 
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A 


Expressing * and @ in terms of Oand de 


q a ee (I-38) 
Fi ©@ im 


After integrating (38) the following results 5 


oe Sel (cos 6, Cos 9s) 


Po att a 


For an infinite vortex 
as 


Wp = ath S20) 


and for a semi-infinite vortex, 


ee i a 
ip a Ghh (I-h1) 


In the present experimentation the value of induced flow at a point P 


will always be in the region 


j7 es fa 
oe s Wp ae (1-42) 
nara ATT A 
ferme une Static case. 
CC 


Since © a > eee 


Then Ces a — a 


Hence from (25) 


Sp Ge (1+ Ces On) 


ATi In 


Now assuming a constant spanwise downwash distribution in the work of 
Spreiter and Sacks (Reference 26), we know that the trailing vortex 
core will flow in the direction of the flap deflection with a velocity 


of approximately 
Aa = at tan Cm e WS 


or as in Reference 18 


W eee caer FI Tan [Se Be, ie \ od VV ><? 
LTR PAL, 


Thus , 


. 
Veore (41) CL 


eon a 
V5? 


TPR A TO 


Then by substituting (T-29) and (I-35a) into (I-46) we arrive at 
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Vena _ tl Spf St GY 
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(I-43) 


(r-44) 
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Transposing (1-47) we arrive at 


Veort a 4] S¢ )¢ V5 4 C V2 +6; | 





a 
ie = - se (1-48) 
7 TPR Yoo” + AA Wy’ 
e 
Assuming 
2 a 
a oe - \f 
Ve > Lh V 
7 FX J? = ie 
and : 
og oe V/ 
a 2 
C, lo >? a (I-50) 
Veore 
feeean readily see that the value of Cer ¢ 


is independent of velocity and is dependent only on vf. This theoretical 
analysis corroborates the previous hypothesis that ohas no affect on 


either or cosé. 


ae [¢ + 1G] a ce 
. Ved Saal 4 
Wie 


Thus, the partial of circulation with respect to '» approaches the value 


given as follows: 


oes 4 ve (I-51) 


For velocities decreasing from 100 ft/sec to 50 ft/sec 


em — e 


- 
Vor (1-52) 
Hence, the partial derivative is a positive number, or 
| ee (155) 
A typical plot is shown in Figure 38, for oe a, Muncy 1on of Lo. The 
actual value, though positive, has a small value and hence the affect 


of oon |! is very small. 


Previously, it has been proven that 


| See | 
oy as : oo = ¢ 
dh a =O (1-54) 
Q Veo 3 VO 
e ° 4 ° Vaoee i. J y ’ 
fi-emcoO 1S Obvious from the plot of ———- as a function of ua? 
Ve 


in Figure 38. 
From (I-43) it was shown that 
— ; aa ae 
Wp P| Pies 2 ay) (I-55) 


Assuming (I-53) is true and substituting (I-54) into (I-55) it is 


obvious that We is a function of V0 but is influenced very weakly by 


a change in hea. Bane lis & 


W 


and is approaching a constant at high yy. 
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Neaaseore 2ntl ence BS on“) can be investigated. 


we know that 


Sy = FLA, cos a 


Readily, we can see from that (I-35) that 


we 2 


i 


BPS - 


constant 
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Geometrically, the value of + can be described as 





b= 


meen Figure 39, 
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Now, 
W= ¥ - D sin oN ier 
eee as ae ae : 
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Again from (I-55) 


(I-56) 


(I-57) 


(I-58) 


(I-59) 


(I-60) 


(I-61 ) 


i: 


Now the cos 2, is directly related to h ia Ghee 





cos o): = (I-62 ) 
\D’* +h 
mae Los olf cope, (1-63) 
then from (1-62), cos 0 = Dfeow CX { COLE 
J ee 
' co eee (1-64) 
\(%s Xeon) + 7 
or 
D/,. 3 of Bie CURE 


Coy Oo, = 


ck: cope tD oe (I- 65) 
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Assuming small angles throughout 
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Hence from (I-65) 
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(ee GS), — peer ae ee en eS a oe ee ee ee a (I-67) 
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Also from (I-61) 
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hz {47 - 7 Decree ee. (I-68) 


Now recalling from (1-6) that 


—ew So 
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s j 
oe ae = 4 I 
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We know that the derivative of | eae respect to dy: is 


(I-69) 


} 
mus, J-has no effect ont, Now recalling (I-29) it is known that 


with all other parameters remaining constant 


Dues 
aa Cén STANT 
J Sp A 


memee, trom chain rule differentiation, the results are, 


O <i cove oe 2 OS, cor€ 2 ( 
————, niin ic oe eee ener 


an ee 


P) OF 7 Ge o ope 


eS ae CONSTANT. 


So by substituting (I-71). and (I-73) into (I-72) the result is obvious. 


3 ot LcoRé¢é 
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ee eo OFA) 
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Hence, 2 cep is linearly proportional to oe 


From both (I-67) and (1-68) 
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(I-71) 


(I-72) 


(I-73) 


(1-74) 


(I-75) 


(I-76) 


(I-77) 
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Both values in (I-76) and (I-77) are not constant. 


Incorporating (I-57), (I-76) and (I-77) into (I-43) we realize 


) Ss ast 
a ee 


that 


(I-78) 


This is contributed to by all three variables defined in (I-56). 


From the previous study we conclude that for the standard flapped 
peer oi) without blowing augmentation that (5-0, >4 dp Vio) whereas in the 
co 2 
blown flap case it is Com zs \ (3 \) 6 5 
| Fy ee 


~~ 


In the regime studied, n<1. At higher velocities it approaches 
in Cc 5, 
Wo = p Veo foc 778 =) 
In summary it can be said that there are two variables of importance 
oe 

in the present study. They are Sand Vo, All other parameters involved 
in the generation of induced velocities at remote points downstream are 
also functions of Se and Vos, And, from the previous discussion it woulld 4-277) 
that Se is the more important parameter in the final magnitude of generated 


turbulence. 
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THEORETICAL DISCUSSION OF DYNAMIC INTERFERENCE 


As discussed in Reference ah, the forces on a flapped airfoil are 
a function of its frequency of flap deflection and its aerodynamic 
properties. The magnitude of induced velocities varies as in Figure 40, 
a vector diagram following the Theodorsen function in magnitude and phase 


changes. 


IGure 40 
As can be seen, the phase lags and the magnitude decreases With increasing 


reduced frequency, K. In this discussion 


k= Ws © 


ee | (I-79) 
50 
In the experimental work the reduced frequencies vary from 
ye 
Pek. Sie (1-80) 


Obviously, with such high values of K, the magnitude is decreased by 
nearly half, for the unidirectionally oscillated flap. 


Two unsteady effects have a significant bearing on the forces and 
moments on the wing. "Virtual mass" and ‘wake effect" are the component 


parts that vary the forces and moments on the wing/flap with increased Ky 


The virtual mass effect has no influence on the strength of the shed 
vortex. It does, however, increase the load on the motor which oscillates 
the flap. By incorporating a motor with a large output, this problem is 


ei rO Consequence. 


1 


The wake effect is a 3-D effect as opposed to the 2-D virtual 
mass problem. Although it can be shown mathematically, that the 
flow-induced by a flap oscillation through a plus and minus ©, < would 
theoretically cancel itself out at very high es. the sinusoidal 
oscillation was attempted for a first run. The final value would be 
zero as opposed to the case shown in Figure 43, in which the flap is 
oscillated in one direction only; in this ease ,'Jp would have a constant 


value when Ne > o>, 


For the symmetrically oscillating case the magnitude did begin to 
decrease when higher frequencies were encountered. However, it is 
obvious that the theory will be altered when the movement of the 
vortex core is taken into account such that at mua = G7 , the 
actual Wy, was still greater than zero, although only a fraction of 
the static or quasi static value. Figure 41, iliustrates the wake 


effect for symmetrical oscillation. 


The law of Biot-Savart gives the incremental induced velocity 
at "P" in terms of circulation strength, " ; distance from the 
vortex core, h; and the difference in the cosines of the angles to 
the ends of the segment. At high frequencies, both positive and 
negatively signed vortex circulation components affect the induced 
velocity at P. Therefore, the actual magnitude in the positive 
direction is less than that of the static or low frequency case. 


Figure 42 shows the static case in which the 
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circulation is assumed constant to infinity and that there is no 


dynamic cancellation. 


Ce (a a 

ane. 

i ’ a eel 
Figure Spee | 

Eyeoscillatang the flap am one direécuion, the ondy real aifeeu on 


the sinusoidal frequency response is to input a bias such that as (Jc SO 


U 


the US» value apuredenes 42ero..) 1ence. assuming at present that the 
vortex core position is constant, the oscillatory component of the 
induced velocity will symmetrically decrease about the mean value of 
, aS VY wD 
(I-79) 


Uncen tell at Wo, Will actually equal this constant value and the 


oscillatory gust will no longer be present. 


Figure 43 illustrates the typical core motion and accompanying 


Circulation strength for a deflection in a single direction. 
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panee the vortex core does cross the plane of symmetry it is obvious 


that the bias will change and that the value at UJACwill be 


/ (I-80) 
oe ates < Wy cee i) 


Eas fees 


The sensor dynamics had little or no affect on the output data since 
the experimental oscillation bandwidth of the present discussion was far 
below the natural frequency of the’sensor. As described in Appendix IY, 

ay = .f for the sensor. These frequency and damping characteristics 
assure fidelity in the measurement of the induced velocities. For 
experimentation at much higher frequencies (i.e.; (Ye > ho hz) sensors 
with correspondingly higher natural frequency must be substituted for 


the present vanes. 


APPENDIX II 
INITIAL EXPERIMENTAL WORK IN TWO DIMENSIONS 


Experimental trends in Ce generation for the blown flap were 
arrived at through a comparison to an analytic and experimental study 
of a simple airfoil section with known aerodynamic properties. Some 


characteristics of the blown flap model had to be preserved. First, 


the slot position was investigated. From a physical point of view, the 


slot must be directed to produce an increase in boundary layer energy 
in the streamvise direction. This re-energization of the boundary 
layer will tend to maintain attached flow over the flap and 

thus from supercirculation adGgitien, Increase ( . 

Again, as previously analyzed, the vector component of momentum also 
adds to the oe The slot is chosen at the flap hinge point since 
it is there that separation would be likely to occur and it must 
operate bi-directionally. Since,in all cases,the forward SbSe ne ee 
point will always cccur exactly at the leading edge tor 2 een 

and Curd excepG Gor very lance me , leading edge blowing ss 
not necessary. Slot area was designed at a width capable of being 


choked at the required flow rate by the Schramm compressor. 


The chord of the airfoil and flap are important in lift and 
circulation development. To correlate the 2-D experiments with those 
of Hazen and Knowlton (References 10 and 14), the overall chord was 
chosen to be 16 inches over the tunnel velocity range, the Re's 
agree well in both 2-D and 3-D cases in concern. The flap size 
is of great importance since the airfoil will always have an : 
imibial \ev hence, the larger the deflected area, the larger the 
Garculation produced. The ITlap size, however, must also be 
considered from an actuaticn load viewpoint. Thus, the required 
circulation must be weighed against the required input power to 
oscillate the flap at high frequency. A basic .45¢e flap seemed 


a reasonable initial compremise. 
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The tradeoff of deflection angle for inertia of the flap 
was also a basic concern. Although very high flap deflection 
angles could produce high (; , it is also a fact that the blown 
air produces dynamic loads which increase with Se and add to 
the already large inertial loads at high frequencies. Previous 
experimentation by Korbacher (Reference 15) with jet flap blowing 
has lead to a group of optimal parameters for a reas model. 
im his work, Korbacher optimizes the lift and drag. In the 
present study, the drag 1s of little importance since it is 
circulation only that produces shed vortices. Thus, alarge 
a deflection as possible would be the optimun, fctuation 
mechanics again limit the range of operations. 
Flap deflection greater than 60° would be difficult to oscillate 
with the present equipment; hence, Dp TO was taken as the 
maximum with variable connecting arm positions available for +45° 


and +30° , 


2-D Smoke Tunnel Testing 


The 2-D smoke tunnel experimentation was used as a means of 
observing trends in, production due to variations in flap chord, 
Se and. Cu . The plot on Figure 47 gives some idea as to the 
Results OL Ghese Cee Knowlton in Reference 1h describes a, method 
of calculation of 4; by streamline deflection upstream. In the 
present study a simple comparison was made between the jet flapped 


models and the 2-D symmetric airfoil with assumed C & = sl/den. 


oul 


2 


Having made this assumption, the 2-D symmetric airfoil without 


blowing was given increments of “X. Observing the streamline 


deflection at a point K upstream 


eee 





of the leading edge in Figure 44. Thus a graph of & versus As 
and hence, _, versus hs Was arrived at for the perfect 2-D, 
PMiigace ahr On J, 


The assumptions made to this point are: 


Cp eer nate a 
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(4) varying Se 


(fo\* 50 ft/sec 


(5) tunnel velocity 3 


The results of these observations appear on Figure 45. The oe. 
versus ae has a constant slope. Thus, this curve can be used 


to determine C for other configurations at o< = 0 but with 
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FIGURE 45. BLOWING AUGMENTATION POWER SETTING. 


Results of this comparative ore ee appear in Figure 46 for various 
CA and. O ¢. at Cos .25¢c and then .35c. As can ve seen, the overall 
trend illustrates that the various partial derivatives in () production are 
all positive. From these conclusions and knowledge of dynamic inertial 


load difficulties, the following values were chosen: 


(3) CHOKED FLOW 


Vortex generation was attempted at NASA Langley (References 7 and 22) 
in the Low speed wind tunnel facility. High frequency oscillation of 
the entire airfoil seemed mechanically unfeasible. By decreasing the 
size of the oscillating area to just the flap aided with blcwing, the 
large (7 Gs possible without the tremendous inertia loads at high 


— e e e ra ~ : 
frequency-As discussed previously from Figure 46, the values of 3 ey 


| d Se 
me 


Ce 
and are obviously high. With Fe increased to .45c and Oy 


yan 
to a choked condition the entire C. will increase greatly. 
~“ 


CONCLUSIONS OF 2-D SMOKE TUNNEL EXPERIMENTATION 


1) A basic symmetrical airfoil with a .45e flap will give an 
adequate tradeoff between Ce, and the aerodynamic and 


inertial load effects at high ie 


2) Choked flow at the blowing slot will give an even emission 
of high velocity air across the span of the model while 


producing the largest possible AC.. 


ror. 
3)  S;= 1G0 will be the cptimal tradeoff deflection between 


Jx{,. and mechanical driving problems. 


4) From the graphs (Figure 44) previously presented of 
experimental e pPreduciiOn, parvial derivacaves Can be 


assumed as follows: 
Ot ‘ ‘) y/ 
— as kh fr 
3 be 


and 


Me o 
© K vf 4 
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APPENDIX III 
SCALING 


In wind tunnel testing, scaling of parameters is very important 
in the simulation of real atmospheric situations. Of major importance 
in producing acceptable turbulence in a 3-D tunnel is the frequency 
of the sinusoidal gusts being generated with respect to the tunnel speed. 
Reduced frequency, K = ior is a nondimensional parameter that 
~ Vos 
expresses this relationship and represents a grouping of the important 


dimensional characteristics. 
To begin,we first hypothesize a full seale aircrart and a model 


Size. If we assume a gust wave length equal to the overall length of 


the aircraft we would have typically 


Pe 
a O21 (Tit 


and 


) aan Gane) 
ie v~ 1 \ 


Next the range of K possible for the wind tunnel can be calculated. 


i 
ra ae yi {> 
K wCoD = @ atta, Nanay A 5 oe (III-3) 
ee 
es 
Here 
n he « 
yy yn.OD en ) 4G = BI ae ) (III-4) 
with } 
ey ee x) Or) Gee, 
a: a vy See <x vo ae. < t oe, J yi a = . Cri t= 5.) 
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and. 


Ou <x i) MOD S | 2 (III-6) 


the possible range of K is 


© Salk = 4 OL base 


Assuming 
Vak = GOO Ft/se ae 
and 
Lal. =A ole = So +t. (IIT-9) 
and since 


on (III-10) 


for a correct simulation, 


eee V | 
Ware >Krop "Ale K 14) hs (III-12) 





or 


co Cor ae ee Z > 
a ae (ieee) 
The high frequency simulation is of interest for the study of 
structural and aeroelastic problems, whereas the low frequencies may be 
of interest in simulating pilot-control inter-relationships in the 


stability and control region of aerodynamic study. 
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APPENDIX IV 


VANE DYNAMICS 


i 








The dynamics of the vanes are very important to accurate 
measurements in both magnitude and frequency. Graham (Reference 8) 
undertook to hypothesize the dynamics of just such a vane with. 


optimal dynamic characteristics. 

Vanes, connecuing arm, and potentiometers are the essential parts 
to the angle measuring system. All contribute to the overall moment~ — 
of inertia of the system and thereby affect the natural frequency 
and damping in the simulated second order system. 

Such a system can be described with the basic equation of motion, 


a a8 a . ae 
to Sy, i ge 1 on, = a Neon CG = C) (TV-1) 
a 


This equation describes the summation of torques due to angular 


acceleration, angular velocity and angular position. 


First, the torque tending to bring the vane to zero angle of attack 


is 
(Iv-2) 


This equation is the summation of the products of lift and moment arms 


fiomesouh tere and afb Lifting surfaces. 


For small angles 


Coe | ae 
ye. © 5) (IV-2a ) 
Bence, : 
a a | | 
eee Th, +Da el Le - hic +Dr 6] Up (IV-3) 


The surface area of both vanes is identical so 


+= Pineda O] | Le Ar | (rv-h) 
Now 
L Cag dd (Iv-5) 
p= \lo*4 Cre | 4 5 (IVv-6) 
Bes cence 


x = © (IV-6a.) 


Thereby 
ie ~ e a 
Se (Iv-7) 
D = Kore ° Ce | 4 _* Gnee)) 
Substituting in the torque equation 
: “i c. r vb a <~ a 4 tes / e : 
wl 8] qv FL Os 0 o8]95s5 | Le de | 
or a ; 
es z 4 cae e oe } / 
O. = © ; } ap (> : ( & | v v 4 Fi “| ; 
{ | = by Ss t DoW_ ~ (Iv-9) 
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Differentiating with respect to O yields 


at = tm OMe, Gal 4 Sy [hi Ae (Iv-10) 


For small angles of attack 


and Cy, is approximated by Cre 


Thu lige ha cia = 
; Lae a | c moa i (y min | q Sy Ae ke | (IV~11) 
Dae, \ 


Next, the torques proportional to angular velocity are important. 
As angular velocity occurs a damping opposing the motion is 


developed from the change in Cee) eee to. the free 


stream. _ = 26 6 


uc 
- ca ) Sc oy ; 9 / 
G2 = |G : iy 25416 tall 4 Tae 
: Vx ae _ 
y a Ap 3 ai od r as Iq Sy YP (lee) 
2 Vise) Uae as 


Differentiating with respect to c we see 


_—s 


(IV-13) 
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Collecting terms yields 
= =(Cy+ min > S \La he | 
aC ve 


(Iv~-14 ) 


Substituting (Iv-14) and (IV-11) into (IV-1) we arrive at the 


characteristic equation of motion for this second order, mass-spring- 


damper type system. 


Té + ue iS + At @ =O | (Iv-15) 


or 
a 7 


Té 3 '@ mls Sy, Tle /6 +k “ it Onin Dy Cetee “rae 


By comparing coefficients of (IV-16) and (IV-1) we see that 


. \¢ he +, min_| Qoy. \ J A -f z{ Givqi) 


and 
| 


For the experimental model being used and the dimensions of the vanes, 


the following dimensional characteristics apply: 
= Bs ~6 | aia, ine 4 
OP 6) Gis, wah 
i V 


€ te ey . i} 
Oe el ae 
‘ 
CMe 
4 ae / 


VALI SL io 
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Solving equation (IV-17) ee 


Um 239 Hz 


and 


The following vane system is therefore critically damped with a natural 


frequency of 35 hz as indicated. 


Experimental tests were made and the vanes react exactly as the 
theory would predict. Figure 49 is an example of the vane reaction 


to a step input. 


: \ : 
It can be concluded that since the '“> of the vanes is far above 
that of the actual testing frequencies and since the vanes are 
critically damped, no overshoot will occur even at high frequencies 


and measurement accuracy in the dynamic case is commensurate with 


that in the static case. 
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